NASA Contractor Keport 4044 


TRANDESNF: A Computer Program 

for Transonic Airfoil Design 

and Analysis in Nonuniform Flow 


J. F. Chang and C. Edward Lan 


GRANT NAG 1-308 
JANUARY 1987 


NASA 


NASA Contractor Report 4044 


TRANDESNF: A Computer Program 

for Transonic Airfoil Design 

and Analysis in Nonuniform Flow 


J. F. Chang and C. Edward Lan 

The University of Kansas Center for Research, Inc. 

Lawrence, Kansas 


Prepared for 

Langley Research Center 

under Grant NAG1-308 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1987 



TABLE OF CONTENTS 


Section Page 

Introduction 1 

Summary of the Theoretical Method 2 

Program Capabilities 3 

Input Data Format 4 

Output Variables 14 

Program Execution . 16 

Sample Input and Output 18 

References 66 


PRECUHNQ PApE H^ANK NOT PILMEf) 



INTRODUCTION 


In this report, the usage of a transonic airfoil code for 
analysis, inverse design, and direct optimization of an airfoil 
immersed in propfan slipstream is described. For a detailed 
description of the theory, Reference 1 should be consulted. 

In the following, a summary of the theoretical method, program 
capabilities, input format, output variables, and program execution 
are described. Input data of sample test cases and the 
corresponding output are given. 



SUMMARY OF THE THEORETICAL METHOD 


In the present method, the Euler equation is simplified by 
introducing a velocity function $ and a rotation function F. The 
latter is to account for the effect of flow nonuniformity, while the 
former is similar to the total velocity potential in the potential 
flow theory. The resulting equation can be shown to be (Ref. 1) 

(a^ - u^)<& - 2uv$ + (a^ - v^)$ = uvF - (a^ - u^)F (1) 

v 7 xx xy yy y x 

where a is the local speed of sound and u, v are total velocity 

components in the x and y directions, respectively. Since the left- 

hand side of Equation (1) is the same as the full potential 

equation, it is convenient to modify an existing full-potential 

transonic airfoil code to solve the present problem. For this 

purpose, Carlson's code (Ref. 2) was chosen. The method of 

Reference 2 is based on a finite-difference approximation of 

Equation (1) in a Cartesian coordinate system. The finite- 

difference equations are solved by column relaxation. In the 

present problem, the rotation function and the stream function are 

also determined in the relaxation process. 

To design an airfoil by direct optimization, the airfoil shape 

is expressed in a Fourier cosine series internally in the code with 

the Fourier coefficients being the design variables. Optimization 

of lift-to-drag ratio is accomplished by using CONMIN optimizer 

(Ref. 3) with lift and trailing-edge closure constraints. 
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PROGRAM CAPABILITIES 


This program has the following features: 

(1) It is applicable to the analysis, inverse design, and direct 
optimization of an airfoil in a transonic uniform or nonuniform 
flow. Boundary layer calculation may be included. 

(2) The nonuniformity may be prescribed in the form of freestream 
Mach profile and/or temperature profile. 

(3) The airfoil may be located vertically at any place in the 
nonuniform region. 
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INPUT DATA FORMAT 


*** ALL INPUT DATA ARE IN THE LIST-DIRECTED FORMAT. *** 

GROUP 1: CASE DESCRIPTION 

TITLE: DESCRIPTION OF THE RUN CASE 

GROUP 2: READ INPUT OPTION 

IOPT EXECUTION OPTION: 

= 1 ANALYSIS OR INVERSE DESIGN 
= 2 PLOT LINEAR, NOT AVAILABLE IN THIS VERSION 
= 3 PLOT CONTOUR, NOT AVAILABLE IN THIS VERSION 
= 4 OPTIMIZE 

= 5 OPTIMIZE WITH PLOTS, NOT AVAILABLE IN THIS VERSION 


IF IOPT = 1, SKIP GROUPS 3-15. 

GROUPS 3-15 ARE CONMIN PARAMETERS (SEE REF. 3). 


GROUP 3: 

IPRINT = 0 PRINT NOTHING 

= 1 PRINT INITIAL AND FINAL INFORMATION 
= 2 FIRST DEBUG LEVEL 
= 3 SECOND DEBUG LEVEL 
= 4 COMPLETE DEBUG 

ITMAX MAXIMUM NUMBER OF ITERATIONS IN THE OPTIMIZATION 
PROCESS 

NSCAL SCALING CONTROL PARAMETER 

= -1 IF USER SUPPLIES SCALING VECTORS. 

= 0 IF THERE IS NO SCALING. 

= +1 FOR AUTOMATIC LINEAR SCALING EVERY NSCAL ITERATION. 

GROUP 4: 

NFDG GRADIENT CALCULATION CONTROL PARAMETER 

= 0 IF ALL GRADIENT INFORMATION WILL BE CALCULATED BY 
FINITE DIFFERENCE. 

= 1 IF ALL GRADIENT INFORMATION IS PROVIDED BY SUBROUTINE 
ANALYSIS. 

= 2 IF GRADIENT OF OBJECTIVE FUNCTION IS PROVIDED BY 
ANALYSIS. 

USE 0 IN THIS VERSION. 
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GROUP 5 


FDCH = 0.01 TO 0.025 TYPICALLY FOR TRANSONIC AIRFOIL DESIGN. 
RELATIVE CHANGE IN DECISION VARIABLE IN CALCULATING 
FINITE DIFFERENCE GRADIENTS. 

FDCHM = 0.00075 TO 0.0005 TYPICALLY FOR TRANSONIC AIRFOIL 

DESIGN. MINIMUM ABSOLUTE STEP IN FINITE DIFFERENCE 
GRADIENT CALCULATIONS. USE A SMALLER VALUE WHEN THE 
AIRFOIL IS CLOSE TO THE OPTIMAL SHAPE. 


GROUP 6: 


ALPHAX = .1 TYPICALLY 
ABOGJ1 = .1 TYPICALLY 

GROUP 7: 


NCON = NUMBER OF CONSTRAINT FUNCTIONS. 

NSIDE = SIDE CONSTRAINT PARAMETER 

= 0 IF DESIGN VARIABLES DO NOT HAVE LOWER OR UPPER 
BOUNDS. THIS IS THE OPTION IN THIS VERSION. 

= 1 OTHERWISE. 


GROUP 8: 


LINOBJ = LINEAR OBJECTIVE FUNCTION IDENTIFIER. 
= 0 FOR NONLINEAR OBJECTIVE FUNCTION. 

= 1 FOR LINEAR OBJECTIVE FUNCTION. 


GROUP 9: 

CT = -.1 TYPICALLY. CONSTRAINT THICKNESS PARAMETER. 

CTMIN = .004 TYPICALLY. MINIMUM ABSOLUTE VALUE OF CT 
CONSIDERED IN THE OPTIMIZATION PROCESS. 

CTL = -.01 TYPICALLY. CONSTRAINT THICKNESS PARAMETER FOR 
LINEAR AND SIDE CONSTRAINTS. 

CTLMIN = .001 TYPICALLY. MINIMUM ABSOLUTE VALUE OF CTL 
CONSIDERED IN THE OPTIMIZATION PROCESS. 

THETA = 1. TYPICALLY. MEAN VALUE OF THE PUSH-OFF FACTOR IN 
THE METHOD OF FEASIBLE DIRECTIONS. 

PHI = 5. TYPICALLY. PARTICIPATION COEFFICIENT, USED IF A 
DESIGN IS INFEASIBLE. 


GROUP 10: 


ITRM = 3 TYPICALLY. NUMBER OF CONSECUTIVE ITERATIONS TO 

INDICATE CONVERGENCE BY RELATIVE OR ABSOLUTE CHANGES, 
DELFUN OR DABFUN. 
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GROUP 11 


DELFUN = .001 TYPICALLY. MINIMUM RELATIVE CHANGE IN THE 
OBJECTIVE FUNCTION TO INDICATE CONVERGENCE. 
DABFUN = .001 TYPICALLY. MINIMUM ABSOLUTE CHANGE IN THE 
OBJECTIVE FUNCTION TO INDICATE CONVERGENCE. 


GROUP 12: 


N1 = NDV + 2, WHERE NDV = NCOEF*2 

NCOEF NUMBER OF TERMS IN THE FOURIER SERIES TO REPRESENT THE 
AIRFOIL UPPER AND LOWER SURFACES. YJ TYPICALLY. 

IREP = 1 FOR READING DESIGN VARIABLES FROM PREVIOUS RUN FILE. 
RECOMMENDED IN SUBSEQUENT ADDITIONAL RUNS. 

= 0 OTHERWISE. 


IF NSCAL £ -1, SKIP GROUP 13. 


GROUP 13: 

SCAL(N5) VECTOR OF SCALING PARAMETERS, 40 VALUES. 

TYPICALLY, 40*0.001 FOR AIRFOIL DESIGN. 

GROUP 14: 

UPLIFT UPPER LIFT CONSTRAINT 

DNLIFT LOWER LIFT CONSTRAINT 

THTE TRAILING-EDGE THICKNESS 

GROUP 15: 

XI THE X-COORDINATE BEYOND WHICH THE AIRFOIL SHAPE IS TO 
BE OPTIMIZED. NOTE THAT X(LE) = -0.5 AND X(TE) = 0.5. 


*** THE FOLLOWING ARE MOSTLY TRANDES PARAMETERS. SEE NASA CR-2821 
BY L. A. CARLSON. *** 


GROUP 16: FREE STREAM CONDITIONS 


M MACH NUMBER 

ALP ANGLE OF ATTACK IN DEGREES 

CIR NONDIMENSIONAL CIRCULATION (HALF LIFT COEFFICIENT); 

USUALLY STARTS WITH 0.0 

CDCORR WAVE DRAG CORRECTION AT ZERO ANGLE OF ATTACK. 

CORRECTION SHOULD BE DETERMINED FOR EACH AIRFOIL 
AND GRID COMBINATION. 

RN REYNOLDS NUMBER 
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GROUP 17: COORDINATE STRETCHING PARAMETERS. 

DO NOT CHANGE THE SUGGESTED VALUES OF THIS GROUP. 

A1 STRETCHING CONSTANT FOR THE Y DIRECTION. SUGGESTED 
VALUE 0.246. 

A2 X DIRECTION STRETCHING CONSTANT IN THE AIRFOIL REGION. 
SUGGESTED VALUE = 0.15. 

A3 X DIRECTION STRETCHING CONSTANT OUTSIDE THE AIRFOIL 
REGION. SUGGESTED VALUE = 3.87. 

A4 THE POSITIVE X LOCATION WHERE THE COORDINATE STRETCHING 
CHANGES. SUGGESTED VALUE = 0.49. 

S4 CORRESPONDING X4 VALUE IN THE TRANSFORMED PLANE. 

SUGGESTED VALUE = 2.0. 

GROUP 18: OPTION SPECIFIERS 

IUNIFM FLOW TYPE INDICATOR 

= 0 FOR NONUNIFORM FLOW 

1 FOR UNIFORM FLOW 

INV PROGRAM MODE PARAMETER 

= 0 FOR ANALYSIS OR DIRECT OPTIMIZATION 

1 FOR DESIGN 

ITACT VISCOUS INTERACTION INDICATOR 

= 0 FOR INVISCID ANALYSIS AND DESIGN 

1 FOR ANALYSIS WITH INTERACTION 

ITERP INTERPOLATION INDICATOR IN INVERSE DESIGN FOR 
INTERPOLATING INPUT CP DISTRIBUTION IN GRID 4 
= 0 CP TO BE READ IN GRID 4 
= 1 CP INTERPOLATED FROM GRID-3 INPUT 

IDEBUG DEBUG INDICATOR 
= 0 NO DEBUGGING 

1 DEBUGGING 

2 CHECK PRESSURE DISTRIBUTION AT THE STABLE 
REGION FOR DESIGN 

IREAD INDICATOR TO USE RESTART DATA 

= 0 IF RESTART DATA ARE NOT USED 

1 IF RESTART DATA ARE USED AND CONVERGED SOLUTION 
WILL BE COPIED BACK TO DISK. 

= 2 IF RESTART DATA ARE READ FROM, BUT NOT WRITTEN ON, 
DISK 32. IN DIRECT OPTIMIZATION, DISK COPYING 
WILL BE DONE AT THE END OF EXECUTION. 

I TIME PROGRAM RESTART INDICATOR 

= 0 FOR NORMAL RUN 

1 FOR 10-MINUTE CPU LIMIT BUILT IN THE PROGRAM 

I HALF OPTION TO DETERMINE STARTING PHI VALUES FOR THE FINE 

GRIDS FROM VALUES CALCULATED IN THE COARSE GRIDS. 

= 0 AVERAGE VALUE 

1 ASYMMETRICAL QUADRATIC INTERPOLATION 

2 ASYMMETRICAL CUBIC INTERPOLATION 

3 SYMMETRICAL CUBIC INTERPOLATION 

NOTE: IF IUNIFM = 1, SET ITERP = 0. TYPICALLY, I HALF = 0. 
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GROUP 19: 


NUMBER OF SPECIFIERS 


IMAX NUMBER OF VERTICAL GRID LINES IN THE HORIZONTAL 
DIRECTION. LIMIT TO 97. 

JMAX NUMBER OF HORIZONTAL GRID LINES IN THE VERTICAL 
DIRECTION. LIMIT TO 97. 

NHALF NUMBER OF GRID REFINEMENTS. TYPICALLY, 2 
MITER MAXIMUM NUMBER OF ITERATIONS. TYPICALLY, 1600-3200 
LP RELAXATION CYCLE INTERVAL IN WHICH DETAILS OF BOUNDARY 

LAYER SURFACE COORDINATES, ETC., ARE PRINTED OUT. 
TYPICALLY, 4000-6000. 

GROUP 20: VISCOUS INTERACTION PARAMETERS 

INPUT TYPICAL VALUES FOR NO VISCOUS INTERACTION CASES 

ITEUPC UPPER SURFACE TRAILING EDGE CORRECTION INDICATOR. 

ONLY USED IN THE VISCOUS INTERACTION CASE. 

= 0 FOR NO CORRECTION (TYPICAL) 

1 WITH CORRECTION 

ITELWC LOWER SURFACE TRAILING EDGE CORRECTION INDICATOR. 

SAME AS ITEUPC. 

SP MAXIMUM VALUE ALLOWED FOR THE NASH-MACDONALD 

SEPARATION PARAMETER WHEN X < XSEP. 

TYPICALLY, 0.004. 

XSEP X LOCATION AFTER WHICH NASH-MACDONALD SEPARATION 
PARAMETER CAN ASSUME ITS CALCULATED VALUE. 

TYPICALLY, 0.44. 

XLSEP LOCATION AT WHICH THE TRAILING EDGE CORRECTION 
PROCEDURE BEGINS. IT SHOULD CORRESPOND TO THE 
POINT OF SEPARATION. TYPICALLY, 0.50. USED ONLY IF 
ITEUPC AND/OR ITELWC = 1. 

XPC LOCATION AFTER WHICH LOWER SURFACE DISPLACEMENT 

THICKNESS IS REQUIRED TO CONTINUE DECREASING ONCE IT 
HAS STARTED TO DECREASE. 

= 0.1 FOR AFT-CAMBERED AIRFOILS 
0.5 FOR CONVENTIONAL AIRFOILS 
XIBDLY THE X-LOCATION AT WHICH TRANSITION IS ASSUMED TO 
OCCUR. TYPICALLY, -.44. 

GROUP 21: NONUNIFORM FLOW PARAMETERS 

ISWIRL = 0 NO SWIRL 

1 SWIRL EXISTS DUE TO PROPELLER 
IANALY = 0 UNIFORM FLOW 

1 MACH NUMBER PROFILE BY EXPONENTIAL EXPRESSIONS 

2 MACH NUMBER PROFILE FROM INPUT FOR NONUNIFORM FLOW. 
ISTAG = 0 UNIFORM STAGNATION TEMPERATURE PROFILE 

1 UNIFORM AMBIENT TEMPERATURE PROFILE (TYPICAL) 

2 INPUT STAGNATION TEMPERATURE PROFILE 

3 EXPONENTIAL STAGNATION TEMPERATURE PROFILE 
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IBU = 0 LIFT COEFFICIENT BASED ON FREE-STREAM DYNAMIC 

PRESSURE 

1 LIFT COEFFICIENT BASED ON MAXIMUM FAR FIELD DYNAMIC 
PRESSURE 

2 LIFT COEFFICIENT BASED ON FAR FIELD DYNAMIC PRESSURE 
CORRESPONDING TO AIRFOIL STREAMLINE 

GROUP 22: CONVERGENCE FACTORS 


ICONV CONVERGENCE OPTION INDICATOR. A TYPICAL CHOICE IS 3. 

= 0 CONSTANT CONV FOR EACH GRID 

1 TWICE THE INPUT CONV FOR GRID ONE 

2 FIVE TIMES THE INPUT CONV FOR GRID ONE 

3 FIVE TIMES THE INPUT CONV FOR GRID ONE AND 

TWICE THE INPUT CONV FOR GRID TWO 

4 SAME AS 3, BUT FIX AND RESET RELAXATION FACTOR AND 
SUPERSONIC DAMPING FACTOR AFTER 800 AND 1600 
ITERATIONS . 

5 SAME AS 4, BUT INCREASE ICYL TO 15 AFTER 2400 
ITERATIONS 

CONV CONVERGENCE CRITERION. IN DIRECT OPTIMIZATION, CONV IS 
INCREASED INTERNALLY BY 2.5 TIMES TO REDUCE CPU 
TIME IF CONV < 2.0 X 10 -6 
- l.OE-6 (TYPICAL) FOR NONUNIFORM FLOW 

1.0E-4 (TYPICAL) FOR UNIFORM FLOW 

RC REFERENCE CONVERGENCE FACTOR. 

= 1000. TO 5000. TYPICAL FOR ANALYSIS 
200. TYPICAL FOR DESIGN 
FD RELATIVE CONVERGENCE FACTOR FOR NEWF. 

= 10. TO 50. TYPICALLY. USE 1.0 FOR UNIFORM FLOW. 

SD RELATIVE CONVERGENCE FACTOR FOR SHAPE 

=0.5 TYPICALLY. USE 1.0 FOR UNIFORM FLOW. 

VD RELATIVE CONVERGENCE FACTOR FOR VIS ACT 

= 2. TO 5. TYPICALLY 


GROUP 23: 


RELAXATION FACTORS 

INPUT DEFAULT VALUES FOR THOSE WHICH ARE NOT NEEDED 


W RELAXATION FACTOR FOR PHI. 

1.4 TO 1.7 TYPICAL 

UW RELAXATION FACTOR FOR SHAPE. 

0.4 < UW < 1.0 

RDEL RELAXATION FACTOR FOR BOUNDARY LAYER DISPLACEMENT 

Thickness, used only for viscous interaction for 

IMAX LE. 55. TYPICALLY, 0.20. IT IS REDUCED 
AUTOMATICALLY BY HALF IN THE LAST FINE GRIDS. IT IS 
FURTHER REDUCED BY HALF IN DIRECT OPTIMIZATION. 
RDELFN RELAXATION FACTOR FOR BOUNDARY LAYER DISPLACEMENT 
FOR IMAX. GT. 55. TYPICALLY, 0.125. 
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GROUP 24: DAMPING FACTORS 


EPS SUBSONIC DAMPING FACTOR. TYPICALLY, 0.0. 

EPSS STARTING SUPERSONIC DAMPING FACTOR. 

EPSMIN MINIMUM SUPERSONIC DAMPING FACTOR. ROUGHLY EQUAL TO 
SQRT (MAXIMUM LOCAL MACH NUMBER**2 -1.) 

GROUP 25: AIRFOIL CALCULATION REGION 

THIS GROUP SHOULD BE INPUT ONLY ONCE IN ANALYSIS MODE. 
THIS GROUP SHOULD BE INPUT FOR EACH GRID IN DESIGN. 

XI X LOCATION WHERE DIRECT CALCULATION STOPS. 

= 0.5 IN ANALYSIS MODE OR IN DIRECT OPTIMIZATION. 

IN DESIGN MODE ALLOW AT LEAST TWO GRID POINTS AHEAD 
OF IT. -0.38 TYPICALLY. 

X2 END OF INVERSE REGION. 

= 0.5 IN DESIGN MODE 

10000. IN ANALYSIS MODE OR DIRECT OPTIMIZATION. 


IF IANALY SKIP GROUP 26. IN A UNIFORM FLOW, SKIP GROUPS 26-32. 


GROUP 26: EXPONENTIAL VELOCITY PROFILE PARAMETERS. 

QINF = QI * (1. + AC * EXP (1-((Y - YS) / DD)**2)) 

AC DIFFERENCE BETWEEN PEAK NONUNIFORM MACH NUMBER AND FREE- 
STREAM MACH NUMBER. 

YS VERTICAL LOCATION OF THE AIRFOIL RELATIVE TO THE 
NONUNIFORM STREAM CENTER. 

DD : VERTICAL SPREAD OF THE MACH NUMBER NONUNIFORMITY 


IF IANALY 2, SKIP GROUPS 27 AND 28. 


GROUP 27: NUMERICAL MACH NUMBER PROFILE PARAMETERS 

ISTA1 NUMBER OF STATIONS FOR SPECIFYING MACH NUMBER PROFILE 
AC DIFFERENCE BETWEEN PEAK NONUNIFORM MACH NUMBER AND FREE 

STREAM MACH NUMBER. 

GROUP 28: NUMERICAL MACH NUMBER PROFILE. 

YLOC Y LOCATIONS WHERE VELOCITIES ARE TO BE INPUT, ISTA1 

VALUES. NONDIMENSIONALIZED WITH AIRFOIL CHORD LENGTH. 
INPUT IN THE ORDER OF DECREASING VALUES. 

VINF : CORRESPONDING MACH NUMBERS, ISTA1 VALUES. 


IF ISTAG £ 2, SKIP GROUPS 29 AND 30. 
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GROUP 29 


ISTA2 NUMBER OF STATIONS TO SPECIFY STAGNATION TEMPERATURE 
PROFILE. 

GROUP 30: STAGNATION TEMPERATURE PROFILE. 

YLOC Y LOCATIONS WHERE STAGNATION TEMPERATURES ARE TO BE 

INPUT, ISTA2 VALUES. NONDIMENSIONALIZED WITH AIRFOIL 
CHORD LENGTH. INPUT IN THE ORDER OF DECREASING VALUES. 
TMP : CORRESPONDING STAGNATION TEMPERATURE VALUES, ISTA2 
VALUES . 


IF ISTAG t 3, SKIP GROUP 31. 


GROUP 31: EXPONENTIAL TEMPERATURE PROFILE PARAMETERS 

TEM = (1 + TC * EXP(-((Y - TS)/TD)**2)) 

TC DIFFERENCE BETWEEN PEAK TEMPERATURE AND FREE STREAM 
TEMPERATURE 

TS VERTICAL LOCATION OF THE AIRFOIL RELATIVE TO THE 
NONUNIFORM TEMPERATURE PROFILE CENTER 
TD VERTICAL SPREAD OF THE TEMPERATURE NONUNIFORMITY 


IF ISWIRL t 1, SKIP GROUP 32. 


GROUP 32: SWIRL PARAMETERS 

SWANG SWIRL ANGLE IN DEGREES 

SWBL Y COORDINATE BELOW WHICH FREE STREAM CONDITIONS PREVAIL 

SWBU Y COORDINATE ABOVE WHICH FREE STREAM CONDITIONS PREVAIL 

GROUP 33: 


NI NUMBER OF COORDINATE PAIRS TO DESCRIBE THE UPPER SURFACE 
OF THE AIRFOIL. LIMITED TO 110. 


GROUP 34: 

(XI, YI) AIRFOIL UPPER SURFACE COORDINATES (X, Y) 

NONDIMENSIONALIZED WITH CHORD LENGTH. NI PAIRS. 
INPUT FROM THE LEADING EDGE TO THE TRAILING EDGE. 
X(LE) = 0. AND X(TE) = 1.0. 
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GROUP 35 


UPPER SURFACE SLOPES 


DERIX DX/DS OF THE AIRFOIL UPPER SURFACE AT THE LEADING 

EDGE. TYPICALLY, 0.0. 

DERIY DY/DS OF THE AIRFOIL UPPER SURFACE AT THE TRAILING 

EDGE. TYPICALLY, 1.0. 

DERFX THIRD DERIVATIVE OF DX/DS OF THE AIRFOIL UPPER 

SURFACE AT THE TRAILING EDGE. TYPICALLY, 0.0. 
DERFY THIRD DERIVATIVE OF DY/DS OF THE AIRFOIL UPPER 

SURFACE AT THE TRAILING EDGE. TYPICALLY, 0.0. 

GROUP 36: 


NIB NUMBER OF COORDINATE PAIRS TO DESCRIBE THE AIRFOIL 
LOWER SURFACE. LIMITED TO 110. 

GROUP 37: 

(XIB.YIB) AIRFOIL LOWER SURFACE COORDINATES (x, y). 

NONDIMENS IONAL IZED WITH CHORD LENGTH. NIB PAIRS. 
INPUT FROM THE LEADING EDGE TO THE TRAILING EDGE. 
X(LE) = 0. AND X(TE) =1.0. 

GROUP 38: LOWER SURFACE SLOPES. 

DERIXB DX/DS OF THE AIRFOIL LOWER SURFACE AT THE LEADING 
EDGE. TYPICALLY, 0.0. 

DERIYB DY/DS OF THE AIRFOIL LOWER SURFACE AT THE LEADING 
EDGE. TYPICALLY, -1.0. 

DERFXB THIRD DERIVATIVE OF DX/DS OF THE AIRFOIL LOWER SURFACE 
AT THE TRAILING EDGE. TYPICALLY, 0.0. 

DERFYB THIRD DERIVATIVE OF DY/DS OF THE AIRFOIL LOWER SURFACE 
AT THE TRAILING EDGE. TYPICALLY, 0.0. 


SKIP GROUPS 39-40 IF 

1. ANALYSIS CASES OR DIRECT OPTIMIZATION 

2. MHALF = 1 

3. ITERP = 1 AND GRID FINER THAN GRID 3 
THIS GROUP SHOULD BE INPUT FOR GRID 2 AND 3. 

IF ITERP = 0, IT SHOULD ALSO BE INPUT FOR GRID 4. 


GROUP 39: 


CPU SPECIFIED C p DISTRIBUTION ON THE UPPER SURFACE AT GRID 

POINTS. INPUT FROM THE LEADING EDGE TO THE TRAILING 
EDGE. SEE NASA CR-2821. 
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GROUP 40 


CPL SPECIFIED C p DISTRIBUTION ON THE LOWER SURFACE AT GRIP 

POINTS. INPUT FROM THE LEADING EDGE TO THE TRAILING 
EDGE. SEE NASA CR-2821. 
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OUTPUT VARIABLES 


The following output is available on file it 30 for each grid. 


(1) The input profiles for the free-stream Mach number and 
temperature are first printed. 

(2) Heading 


(3) Listing of input data 

(4) Cartesian grid coordinates. Pairs of I,X(I) are first 
printed. These are followed with pairs of J,Y(J). 


(5) Airfoil coordinates 


X 

YU 

YL 

UPPER SLOPE 
LOWER SLOPE 


HORIZONTAL COORDINATE WITH -0.5 BEING THE LEADING 

EDGE AND 0.5 BEING THE TRAILING EDGE 

UPPER SURFACE ORDINATE 

LOWER SURFACE ORDINATE 

UPPER SURFACE SLOPE 

LOWER SURFACE SLOPE 


(6) Iteration history 

ITER ITERATION NUMBER 

CIR CIRCULATION 

DPM MAXIMUM 4> CORRECTION AT THE GRID LOCATION (I,J) 

NSSP NUMBER OF SUPERSONIC POINTS 

DELTA MAXIMUM BOUNDARY LAYER DISPLACEMENT THICKNESS 

EPSS SUPERSONIC DAMPING FACTOR 

W RELAXATION FACTOR FOR $ 

UW RELAXATION FACTOR FOR AIRFOIL SHAPE IN INVERSE 

DESIGN 


(7) Results of boundary layer analysis (for cases with viscous 


interaction only) 

X HORIZONTAL COORDINATE 

YUORIG ORIGINAL UPPER SURFACE ORDINATE 

DU UPPER SURFACE DISPLACEMENT THICKNESS 

SLU SLOPE OF UPPER SURFACE ORDINATE 

YLORIG ORIGINAL LOWER SURFACE ORDINATE 

DL LOWER SURFACE DISPLACEMENT THICKNESS 

SLL SLOPE OF LOWER SURFACE 
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( 8 ) Pressure distribution on airfoil 

X HORIZONTAL COORDINATE 

CPU PRESSURE COEFFICIENT C ON THE UPPER SURFACE 

CPL PRESSURE COEFFICIENT 08 THE LOWER SURFACE 

(9) Airfoil shape with boundary layer displacement thickness 
included 

X HORIZONTAL COORDINATE 

YU UPPER SURFACE ORDINATE 

YL LOWER SURFACE ORDINATE 

SLU SLOPE OF UPPER SURFACE 

SLL SLOPE OF LOWER SURFACE 

(10) Mach number chart 

"I" increases from top to bottom and "J" increases from 
left to right. The actual value of Mach number is the printed 
value divided by 100 . 

(11) Wave drag coefficient and wave drag correction (CDCORR) 

(12) Plot of results 

UPPER SURFACE C c 
LOWER SURFACE C 
UPPER DISPLACEMENT SURFACE 
LOWER DISPLACEMENT SURFACE 
CRITICAL PRESSURE COEFFICIENT, C * 

LIFT COEFFICIENT FROM CIRCULATION 
LIFT COEFFICIENT FROM C INTEGRATION 
TOTAL DRAG COEFFICIENT P 

PITCHING MOMENT COEFFICIENT ABOUT THE LEADING 
EDGE 

SKIN FRICTION DRAG COEFFICIENT 
PITCHING MOMENT COEFFICIENT ABOUT THE QUARTER- 
CHORD POINT 

(13) CPU time 

In direct optimization, various CONMIN variables may be 
printed, depending on printing options. For details. Reference 
3 should be consulted. 


U 

L 

T 

B 

CPSTAR 

CLCIR 

CL 

CD 

CMLE 

CDF 

CMC4 
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PROGRAM EXECUTION 


This code is written in Fortran 77 language. It is operational 
on the Harris-1000 computer at the University of Kansas and the CDC 
Cyber 175 computer system at the NASA Langley Research center . 

The following files used during execution are defined in the 


program: 

FILE VARIABLE FILE NUMBER USAGE 


INPUT 


JOUT 

TOUT 


LOUT 

KOUT 


20 

22 

30 

31 


32 

33 


INPUT DATA 

TO STORE DESIGN VARIABLES FOR 
RESTART IN DIRECT OPTIMIZATION 
OUTPUT 

DATA IN INPUT-DATA FORMAT 
CONTAINING THE FINAL AIRFOIL 
COORDINATES IN DIRECT OPTIMIZATION 
FOR RESTART 

CONTAIN SOLUTION IN BINARY FORMAT 

FOR RESTART 

OUTPUT FOR DEBUGGING 


File numbers may be redefined in BLOCK DATA. The file #22 is used 
in Subroutines OPT and COEFI. 


Execution in analysis and inverse design is straightforward. 

In direct optimization, the following steps are recommended. 

(1) With an assumed starting airfoil, run the code in analysis mode 
with IREAD = 0 (see GROUP 18). The converged solution is 
automatically saved on file LOUT (#32). 

(2) Change the input file for direct optimization. The following 
values for some input variables are recommended: 


ITMAX = 2 (GROUP 3) 
IREP = 0 (GROUP 12) 
IREAD = 1 (GROUP 18) 
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CONV 


= 2 . 5xlO -6 (GROUP 12) 

RDEL = 0.1 (GROUP 23) 

EPSS =3.0 (GROUP 24) 

(3) The final solution is again saved on file LOUT (#32). If the 
output indicates that the objective function is not changed, 
FDCHM (GROUP 5) may be slightly increased. To restart, copy 
file IOUT (#31) to file INPUT (#20), and IREP (GROUP 12) is set 
to 1. ITMAX (GROUP 3) may be increased. Other variables in 
step (2) remain the same. 

(4) For any subsequent restart, step (3) is repeated. 

(5) After a satisfactory airfoil shape is obtained, file IOUT (#31) 
is copied to file INPUT (#20) again. The file INPUT (#20) is 
then changed for analysis only. For this final analysis, the 
following values for some input variables are recommended: 

I RE AD = 1 (GROUP 18) 

CONV = 1.0xl0“ 6 (GROUP 23) 

RDEL = 0.2 (GROUP 23) 

EPSS =2.5 (GROUP 24) 
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SAMPLE INPUT AND OUTPUT 


1. Input Data for Sample Case 1 


1 


AIRFOIL ANALYSIS IN NON-UNIFORM TRANSONIC FLOW 


.8 

.246 

0 

13 

3 

1 

3 

1 .4 
0 . 


0 

13 

0 

2 


0 . 0 . 

.15 3.87 

1 0 0 
2 3200 6000 
.004 .44 

1 0 

.000001 2000 . 


0 . 

.49 
0 0 


20950000 . 

2 . 


.5 


.52 

.1 

-.3 

-.7 
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3 . Input Data for Sample Case 2 


ESIGN OF AN AIRFOIL IN NON-UNIFORM TRANSONIC FLOW 
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BEGINNING EPSS= 1.40003 BEGINNING Rr>EL = 0.05000 

ACCEPTA3LE REGION FOR ANALYSIS 
AT ITERATION = 1 

EEGIN N GW F 
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0.27024 0PM = 0.0000025 IT 24 21 NSS° = 572 0ELTA = 0.0120 E»SS = 1.40 W = 1.70 UW 

WAVE CD = 0.057925 CDC0RR = 0.000000 

PRESSURE COEFFICIENT 




o 



O 



O 




o 



o 



o 




• 



• 



• 




r- 



M 







II 



II 



ii 




3 



3 



3 




O 



O 



=> 


't 


O 

St 


o 

m 


O 

AJ 

N- 


f— 

O' 


N- 

O' 


A- 

in 

Kl 


• 

Al 


• 

Kl 


> 

AJ 

to 


r- 

K> 


•— 

K1 




• 



■ 



• 



■ 

o 


II 

O 


II 

O 


11 

O 



2 



3 



3 

1 



T— 



O 



T" 




St 



Lf 



>4 


N 


• 

II 


• 

II 


• 

II 



T— 



r- 



r* 


*4 



st 



'f 



st 

o 


II 

U 


ll 

O 


II 

o 

X 



JC 



XI 



£ 

u 


in 

u 


in 

o 


CO 

o 



xn 



m 



to 


>o 


a_ 

oo 


CL 

'O 


ft 

s4 

IA 


Ui 

m 


UJ 

A- 


UJ 

00 

or 



-o 



O 



'O 

N 



A* 


M 

A« 


h- 

A- 

o 


T- 

>40 


Ad 

mo 


r- 

r-O 

o 


T“ 

AJO 


T~ 

mo 



OO 

• 


a 

AJ • 


o 

<4 • 


O 

AJ • 

o 


■ 

mo 


• 

mo 


• 

mo 



o 

■ 


o 

• 


o 

• 




o 



o 



o 

II 


ii 

II 


II 

ii 


II 

II 

u. 



u. 


< 

u 


« 

u. 

o 


m 

o 


M 

Ci 


M 

o 

w 


-i 

11 O 


_J 

II o 


-J 

II O 



UJ 



UJ 



UJ 


00 


a 

<XK- 


o 

CL >4 


o 

cr cn 

Kl 



MO 



HS 



Ms4 

fO 


O 

Ur 


O 

om 


AJ 

Ur 

IA 


A- 

-jm 


A- 

Jtn 


N- 

-Jin 

• 


IA 

u ■ 


m 

O 1 


m 

o • 

O 



o 



o 



o 



II 



li 



ii 

1 



QOH 



CLOM 



ft-OM 




l/>OZ 

-4 


WOZ 

m 


coo z 

m 

II 


VOID 

A II 


WOUJ 

e> ii 


COO LU 

A- II 



ZOM 

O 


ZOM 

o 


ZOM 

o 

Ui 


oo 

AJUJ 


Oo 

AJl-l 


Oo 

AJUJ 

-J 


OOh 

• -J 


mOm 

• -J 


S40MI 

• -J 

3T 


■4 #U 

or 


AJ • U 

o». 


AJ ilk 

or 

U 

O 

Ou. 

1 U 

o 

Ou. 

1 <-> 

o 

OU. 

1 o 


o 

a- tu 


o 

r- UJ 


o 

K» UI 



o 

K> o 

r— 

o 

AJ O 

r- 

o 

>4 O 

m 

*) 

m 

ii i -» 

IA 

m 

II o 

m 

in 

II O 

-n 

N 

O 

m 

K. 

o 

M 

ro 

o 

M 

o 

»A 

• 

<a ui 

11 <4 

■ 

«* cr m 

II -o 

■ 

rfcrui 

li m 

«o 

o 

ac a 

>© 

o 

O COL 

o 

o 

CL ft. 

■o 

O 


inOO 

OO 


won 

ft o 


s4 00 

ao 

■ 


(MOt/5 

c • 


Ad o l/l 

< • 


AJ 04/1 

< ■ 

o 


Otnt/i 

MO 


new 

MO 


DOW 

MO 


II 

OUUJ 

in 

II 

OOIU 

CO 

1) 

OOIU 

v> 


_l 

o a 

ft. 

_J 

o tr 

a 

-J 

O or 

CL 


ui 

o a. 

LJ 

UJ 

o a. 

o 

UJ 

O CL 

O 

li 

c» 

Ost 

II 

o 

Or 

ii 

n 

OO 

II 


tt 

»T— 


Of 

• -O 


a: 

• OU 


a 


ON 

fi 


OKI 

o 


Oft' 

o 

u 

e> 

•o 

LJ 

o 

CO 

o 

o 

A- 

o 


Z 

in 


z 

in 


2 

in 


o 

A— A 

O 

AJ 

H 

O 

r~ 

M 

O 

h- 

IA 

Z 

li • 

T— 

z 

II • 


Z 

II • 

ce 

00 

Z 

O 

>o 

z 

O 

O' 

z 

O 

in 

K 

w 

y 

A- 

H 

3-: 

fs. 

M 

r 

fs. 

• 

o 

a. 

• 

e> 

a 

• 

o 

a. 

■ 

O 

UJt/> 

a 

O 

Ul(/1 

o 

o 

UICO 

a 

O 


DDW 

ii 


CUM 

ii 


CCH 

M 



W 

ro 


C/1 

rn 


to 

r- 



>- 

t\IO 


>- 

N- Q 


M 

oo 



O-J 

r-O 


O-J 

AJO 


o_j 

oo 



o< 

-O 


o< 

A- 


oc 

o 


II 

oz 

AJUJ 

II 

oz 

AJLU 

II 

oz 

AJUJ 

II 

O 

o< 

• > 

o 

r-*< 

■ > 

o 

o< 

•> 

o 

—XT) 

si- 

o*<r 

-JO 

sf 

o<r 

-JO 



-JO 

UO 

oaf 

3 

oo 

• cc 

2 

oo 

•cx 

3 

oo 

O 

rOr 

o o o 

O 

rOr 

o o o 

o 

rOr 

O O O 

o 

o 

u 

r - r* r - 

T“ 

U- 

r- r r 

o 

u. 

T“ r r 

r* 

s» 



"4 



sf 



s4 

• 

Z 

II 

* 

Z 

II 

• 

Z 

M 

• 

T- 

II o 


T" 

II o 


T- 

II o 


M 


Wh 

(X 


COM 

a. 


COM 

tx 



too II 

II II II M 


WO l| 

II II II M 


COO II 

II II II M 



<lu> 

z u 


CLUJ 

Z o 


CLUJ 

Z O 


II 

Ulft z 

7HZOZ 

II 

UJOfZ 

ZM2TOZ 

II 

UlftZ 

ZMZOZ 

II 


o 

OuOhON 


o 

OU0H04 


o 

O O O MOO 


IA 

em)HU.W<WOHin 

in 


CO 

OWHL HCHWH-f 

IA 

in 

ZJK3H-MhliJh00 

Vi 

Z_iM3MC0MUJMft' 

CO 

ZJhZhWMUJM^ 

CA 

a 

m m <• u> *r m «r * 

a 

HfDtfUKHta* 

ft 

Mtr«UJ^M«ror<t 

a 

Ui 

Z-tCCZtt >CC Of 

Ui 

z«<azft>a: cx 

UJ 

z<aza>a: tr 

LU 


ZMUi 

UJ UlOiltJ II 


?Ulu 

Mi ItJlJJIH II 


ZMUJ 

HI UJUlUJII 


-J 

mCLm 

ZhZKJh 

-J 

mclh: 

ZhZhJh 

-J 

mCLMZMZMJM 

-J 

< 

OWHHWHMfr,HOr 

<r 

UIUJHWHHMfOHir 


WlDMHHMHfflMOr 


z 

luu 

O U> «C UJ 

z 

luu o O «*C uj 

z 

UtO 

19 o < ui 

z 

V* 

fPUhlbhU'hl-l-h 

M 

tr oMUJMU'MMMM 

1 -l 

CT' OM UJ MU JMM MM 

M 

u. 

«c <. a. < a. in «i m 

U. 

mc < oj <, iu <c m < h 

u. 

<* <cr.<tu<(/)«LH 

U. 

O «■“ 


o 

O i 

r- 

O 

O 1 

r- 

o 

O ' 


w 


Of 


o 

u 

Ui 
UJQC 
X 

M 

MO 1 
CJCE 
<W 

flDKI 
WE O 

O 

Z 

MAJ 

UJZ UJ 

ftTH or 

« 
ft 

UJHt- Ui 

ft OO ft 

UJZ UJ 

y r> y 

hu h 


54 


VIOLATED CONSTRAINTS 


CTO 


oo 

OO 

OO 





mm 






OO 


OO 


OO 


OO 

oo 

oo 





oo 






oo 


OO 


4 4- 


■4 4 

4 4 

44 





i I 






44 


4 4 


UJ Hi 


UJ Ui 

UiUJ 

UiUJ 





UJ UJ 






UJUJ 


UJUJ 


mm 


inn- 

f 4 

mm 





mm 






mm 




•-in 


oorg 

K>in 

oo 





O' CO 






a-m 


4)43 




m*- 

OKI 

KIK* 





K)4> 






>40 


mm 


K»e 


om 

mo 

mm 





CO O' 






m>4 


mm 


00 OO 


too 

Or- 

mm 





mm 






CM CO 


4343 


ION 


t-4 

r-4 

r-r- 





r-m 






mm 


mm 


• * 


a • 

• a 

• a 





• a 






• a 




OO 


00 

1 1 

oo 

OO 

1 





00 

1 






OO 


OO 

1 


OO 


oo 

OO 

OCT 





mm 






OO 


OO 


OO 


oo 

oo 

OO 





oo 






oo 


oo 


4 4 


4- 4- 

4 4 

44 





( i 






4 4 


4 4 


UU 


UiUJ 

UiUJ 

UiUJ 





UiUJ 






uu 


UU 


>0*0 


r-Ki 

om 

mm 





O.K- 






43*0 


r-r- 


eo 


43 m 

rum 

OO 





mm 






40 


434) 


ino 


oe 

-OM 

mm 





mm 






mo 


mm 


O'T- 


om 

r-O 

Kim 





oo 






O'*- 


mm 


-tiT\ 


oo 

OO 

rgrg 





moo 






>4in 


4)43 


T*r 


*— nj 

mm 

rr 





Mr 






a-*— 


mm 


• • 


a • 

a a 

a a 





a • 






• f 


• a 


OO 

1 


oo 

00 

1 1 

oo 

* 





00 

1 






00 

1 


00 

1 


a-O 


r-O 

r-O 

OO 





mm 






a-O 


OO 


OO 


oo 

OO 

oo 





oo 






OO 


oo 


1 4- 


1 4- 

1 4 

44 





i 1 






1 4 


4 4 


UJUJ 


UiUJ 

UiUJ 

UiUJ 





UJ LU 






UU 


UU 


mo 


43 rg 

-0-0 

mm 





OO 






mo 


r-a— 


r->o 


ON- 

O'*— 

OO 





mo 






me 


4)4) 


em 


Nt» 

ON 

mm 





om 



r- 



4M 


mm 


inm 


om 

<MO 

mm 





r-O 



O 



mm 


mm 


mm 


om 

Om 

rum 





mm 



4 



mm 


43 43 


kiki 


r-K> 

*— K) 

r-r- 





t-4 



U 



mm 


mm 


• i 


• • 

a a 

a a 





a a 



o 



a a 


a a 


OO 


00 

1 1 

oo 

OO 

1 





00 

1 1 



o 

o 



00 

1 


OO 

I 


OCT 


oo 

OO 

oo 





Kl*- 



o 

m 



OO 


CTO 


OO 


oo 

oo 

oo 





OO 



a 



OO 


OO 


4- 4 


4 4 

4 4 

4 4 





1 1 



O 



44 


4 4 


Lulu 


UiUJ 

UiUJ 

UiUJ 





Ui u 






UU 


UU 


Kit- 


me 

om 

mm 





rnr- 






on*- 


r-r- 


O'lA 


COe 

-40 

oo 





40 






om 


4343 


4)-4 


r-in 

mm 

mm 





-4 a— 



II 



43 «4 


mm 


OO 


mm 

-4 1> 

mm 





om 






O' O' 


mm 


on/ 


r-4> 

r-4> 

mm 





oo 



M 



om 


4143 


K>r- 


T*t“ 

r-r- 

r-r- 





4a- 



X 



Mr 


mm 


• • 


• • 

a a 

a a 





• a 



O- 



9 9 


9 a 


00 

1 


oo 

00 

1 1 

00 

1 





00 

1 






00 

1 


00 

1 


OO 

to 

a-o 

T-O 

OO 

O 




r-r- 

m 





OO 

in 

CTO 


OO 

M 

oo 

OO 

oo 

O 




OO 

o 





oo 

M 

OO 


4- 4- 

z 

1 4 

1 4 

4 4 

4 

O 



1 1 

i 


m 



4 4 

Z 

4 4 


LULU 

M 

UiUJ 

UiUJ 

UJU 

Ui 

O 



UIUJ 

UI 


o 



UU 

M 

UJ LU 


K1K 

c 

men 

me 

mm 

a— 

1 


cr 

Om 

in 


l 



mm 

C 

r-r- 


K.r- 

a 

00 a- 

me 

OO 

m 

Ui 


o 

mm 

m 


u 



Mr 

O- 

4343 


em 

M 

t> N- 

Om 

mm 

Ui 

on 



m*~ 

to 


4) 



•4 m 

M 

mm 


oo-* 

m 

a— 4) 

40- 

mm 

^»m 

o 


z 

inO 

Kl 


T— 



CO 4 

m 

mm 

r*. 

43 e 

Z 

4)f\i 

mm 

mm 

045 

m 


o 

O O 

m 


■4 



4>-4 

z 

4343 

U 

n^co 

o 

a~ m 

a— eg 

r-r- 

CO 

4) 



Nr 

T” 


43 


tn 

•4m 

o 

mm 

c 

• a 

u 

• a 

a a 

• ■ 

Z a 

m 


Ui 

a • 

a 


•4 


M 

9 9 

4-1 

a a 

z 

OO 


oo 

OO 

OO 

SO 

a 


13 

OO 

O 


a 

in 

Z 

OO 


OO 

• 

1 

o 

1 

| 

1 

r— 

O 

r* 

Z 

1 

1 


O 

M 

M 

1 

o 

1 

a— 


LU 




li 

1 

O 

c 




1 

Z 

C 


u 


H 

nr* 

M 

• r 

*~t- 

OO 

MO 


4 

X 

r-r— 

a- 



M 

a 

VI" 

M 

OO 

M 

OO 

C 

oo 

OO 

oo 

o 


UJ 

o 

aOO 

O 



C 

M 

OO 

c 

00 


1 1 

-J 

1 t 

1 1 

4 4 

S4 

II 

o 


ft 1 I 

r> | 


II 

or 

in 

1 1 

-J 

4 4 

V 

LU UJ 

O 

UJ Ui 

UJUJ 

UiUJ 



o 

o 

O UJ U J 

ft-JJ 

m 


M 

z 

Uiul 

o 

UU 


oe 

M 

k>4» 

om 

min 

HKl 

C 

o 

z 

M-om 

Oo? 


_l 

in 

o 

0-4 

L* 

rr 

M 

•oo 

> 

oo 

mm 

OO 

'-'N- 

M 

o 


Ur O 

MO 


M 

Z 

4-> 

•OO 

> 

4343 

w 

OO- 


mm 

a— N- 

mm 

cm 

UJ 

o 


uoin 

003 


U 

o 


oo 


mrj 

C 

-ora 

o 

>1-0 

r-N 

mm 

M4) 

03 

• 


>ooo 

ujg 



4-1 

o 

4>m 

Q 

mm 

M 

MO 

Zr-mm 

mmm 

mmm 

UJ*— 


o 


i mm 

>m 

ft 


‘■J 

IU 

mo 

ZMec 

UI 

moo 

«c 

mm 

mm 

a * *~ 

Xe 

s 



xmm 

1 K> 

UJ 


ujor 

M 

moo 

c 

mm 

X 

• • 


a a 

a a 

a • 

M a 

or 


m 

w « # 

ii • 

n) 


>c 

C 

• • 


a a 

M 

OO 

Ui 

OO 

OO 

OO 

wo 

LU 


O 

OO 

wo 

il 


►H 

-I 

OO 

Ui 

OO 

w 




1 1 

1 


M 

II 

4 

in t 

1 

rs 

r— 

M4> 

o 


> 

1 



mo: 

a 

ac 

\ 

U 


UJ 

u 

in 

z 

O 

woe 

M 


H« 

N 

*-r- 

MU 

ir*r-- 

Ur-r- 

UiOO 

mm 

y 

c 

m 

Jrr 

‘Mr- 


l 

CU 

> 

T**— 

uunn 

WO 

nOO 

uojOO 

aoo 

uoo 

u-O 

< 

X 

■4 

UiOO 

OO 

z 

u 

CM 


nOO 

oncoo 

OCO 

cn 1 | 

<i i i 

r i i 

XT 4 4 

O 1 

« 

CL 

m 

c I I 

_J 1 

o 

o 

r-V 

o 

on i i 

c»: 

4 4 

04 

OLUUJ 

CTunu 

OUJLU 

OUUi 

MU 

c 

_J 

O 

►H UI LU 

cut 

M 

O 

O 


o uu 

Dujuj 

MU 

m e 

u.2TNO«t 

Zf-f 

?irin 

OKI 

a 

c 

m 

cr moo 

>«n 

M 

m 

2T 


m-4 

u^Ta-a- 

UK) 

U-CJ43 

o 

*3 a- 

m-o 

OO 

cm 



r- 

cuom 

a 

c 

-4 



u-O 43 

O 

4)4) 

Ca- 

ONM 

HO"C 

MK>e 

Mmm 

u-N * 

M 

Q 

9 

>gN 

Ml/i 

or 

m 

M 


omm 

Mmm 

UO 

m** 

W2UW 

200 

ZKiKi 

m 


U 

O 

4>e 

zk 

u 

a 

UIZ 

u 

me 

wzms 

O 

hr'O' 

KHKlN 

M*4or, 

Mmm 

ucn 

M 

M 

I 


WN 

M 

O 

ftTM 

cr 

M r-O 

Mr* 4343 

UO 

ZO'K. 

zck>ki 

<KiK1 

at r-r- 

um 

< 

C 


OK.K) 

CO 

M 

1 

cc 

< 

ZCMn 

zemm 

u-O' 

UJ i * 

uiar « a 

Or 9 9 

a: a a 

o ■ 

or 

_J 


M a • 

a * 



O' 


Ui a • 

uiir a a 

O a 

MOO 

w KCO 

MOO 

MOCi 

1 o 

M 

o 

II 

WOO 

MO 

Z 


uMm 

u 

MOO 

MKOO 

1 O 

O 1 

CUn | | 

m 

m | 

X 

to 

u 


M | 

<n | 

M 

II 

ertn 

ft 

a t 

OW 1 

X 

< 

<2 

z 

z 

V) 

2 



O 

2- 

L3 


UiZ 

u 

c 

cz 

w 

cr 

frn 

r> 

r> 

T> 

r> 

«r 

m 

Ui 

f ’ 

ij i 

M 

m 

jT 

.V 

ft “ 


“3 

U 

oo 

W 

u 

CL 

u 

o 

O 

o 

4-1 

CkJ 

4-1 

Mu 

M 

o 

oo 

CL 



r>r* 

o 

nn 

rs 




n « 






AA 


nn 

n 

*-co 


t-OU 

*— oO 

a— uU 

r— 




r-OO 

a- 





a— U-« 


a-«J 

r- 

ww 


ww 

ww 

uv 

w 




WVJ 

w 





WW 


ww 

w 


55 


CONSTRAINT PARAMETER* BETA = 0.14060E+01 
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FINAL OPTIMIZATION INFORMATION 
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